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Abstract 

We present a model for diffusion in a molecularly crowded environment. The model consists of 
random barriers in percolation network. Random walks in the presence of slowly moving barriers 
show normal diffusion for long times, but anomalous diffusion at intermediate times. The effective 
exponents for square distance versus time usually are below one at these intermediate times, but 
can be also larger than one for high barrier concentrations. Thus we observe sub- as well as 
super-diffusion in a crowded environment. 
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I. INTRODUCTION 



One of the important issues in biology is to understand how diffusion is affected by the 
environment. This understanding is needed to correctly describe the passive intracellular 
transport as tMs p.oce. my .egulate ..portant cellular p.opert.es: signal transduction fl, 
gene transcription kinetics of reactions [3|] and regulation of cell polarization 

The interior of biological cells represents a very dense and crowded environment with 
a specific molecular mobility. Intracellular diffusion is hindered by barriers consisting of 
large molecules sometimes even immobile tethered large molecules, binding and collisional 
interactions. One way of interpreting such a system is to view it as a disordered system. 



In general for random walks 
expects anomalous diffusion jfi , 
is proportional to the time t: 



l|] in media with disordered microscopic substructures one 
where the mean-square displacement (Ar(t)^) no longer 



{Ar{tf) = e (1) 

with Cq>0. IfO<a<l then we call the diffusion sub-diffusive, and if a > 1 super- 
diffusive; normal diffusion has a = 1. In biological systems as well as for models of biological 
system, so far only sub-diffusive behavior has been observed [6] except if the transport is 
facilitated or restricted. For example if the diffusion is directed by a motor protein, non- 
random super-diffusion can be observed but not for free diffusion. However, super- diffusion 
has been observed in a two-dimensional complex plasma Qj and in two-dimensional Yukawa 



liquids (10| • Thus super- diffusion does exist and its existence in the cell needs to be discussed. 

Sub-diffusion occurs if the mobility of diffusing particle or molecule is impaired by obsta- 
cles (mobile and immobile) or attractive forces. Under these premises it seems unlikely to 
find super-diffusion. 

From a physical point of view the disorder in the cell has a characteristic time scale. If 
the diffusing particle has diffused long enough, such that the ti me scale has been explored, 
one expects normal diffusion. If there is no characteri stic time (in fractal media) then 
the diffusion is always anomalous [8]. Howev er, this picture does not take into account 
the (im)mobility, collision, attractive forces etc. To study the effect produced by immobile 
as well as mobile barriers we have deve loped a model mimicking the molecularly crowded 
environment with mobile as well as immobile barriers taking into account particular length 
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Tried sites after one million time steps 
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FIG. 1: Illustration in small square lattice: Tried (part a) and visited (part b) sites after one 
million time steps. 



scales. 



II. THE MODEL 



In percolation theory 



12], each site on a large lattice is randomly either occupied, with 



probability p, or empty, with probability 1 — p. For percolative diffusion, a random walk 
starts on an occupied site and then on each time step selects randomly a direction to move 
to. It actually moves one unit distance (lattice constant) into this direction if that neighbor 
site is occupied. Empty sites are prohibited for the walker. For p < pc the walk can extend 
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d log<r*r>/d log t; d = 2, a = 0.5, pc-p = 0, 0.1 , ...0.5 
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FIG. 2: Effective exponent for 8000 walks on 7001 x 7001, a = 1/2, A = 0.01. The concentration 
p of allowed sites decreases from top to bottom: pc — p = (+), 0.1 (x), 0.2 (stars), 0.3 (open 
squares), 0.4 (full squares, 0.5 (circles). 



only over the finite cluster in which it started, while for p > pc ii can diffuse towards infinity 
if it started on the infinite cluster. Right at p = Pc anomalous diffusion takes place with a 
mean square distance increasing to infinity but with an exponent a < 1 [12i] . 

In biological applications, the prohibited sites may be more or less mobile biomolecules. 
Their effect can be taken into account approximately by assuming that also a prohibited site 
i allows the walker to move through, with some probability g^. The reciprocal probability 
1/qi then can be interpreted as the lifetime of the barrier, in the sense that about once 
during that lifetime the barrier moves away for one time step before returning to that site. 
Thus we have still a quenched disorder; with annealed disorder where all lifetimes are the 
same, we have normal diffusion, squared distance proportional to time, with a diffusivity 
reduced by the (slowly) moving barriers. We now assume that the probability distribution 
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d log<r*r>/d log t; d=3, a = 0.5, pc-p = 0, 0.1 , ...0.3 
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d log<r*r>/d log t; exponent 0.5, p=0; two dimensions 
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FIG. 3: Effective exponent for 8000 walks on 401 x 401 x 401 (part a, symbols as in Fig.2) and 
7001 x 7001, (part b); a = 1/2, ^ = 0.01. 5 



30,000 (+) and 8,000 (x) walks, d = 3, p-pc = 0.3 
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FIG. 4: As Figure [3] but A = 0.0001 and better statistics or longer times. 



function /(g) for the is a power law, 



/(g) oc l/g» 



(2) 



with some exponent a between zero and infinity. More quantitatively, for each prohibited 
site we determine, when it is visited for the first time by the walker, a random number r, 
homogeneously distributed between and 1, and then fix for that site i as 



(3) 



with some free parameter A. 



III. RESULTS 



Figure d] illustrates for two dimensions at p = Pc — 0.5 = 0.0927 and exponent a = —1/2 
the results of one walk after one million time steps. Part a shows the set of sites which 
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have been tried at least once, and part b shows those sites which have actually been visited 
inspite of the barriers. After 8 million steps, all sites were tried, and after 64 million steps, 
all sites were visited. One can get anywhere, provided one has enough time. 

For 7001 X 70001 square lattices, where Pc — 0.593, our simulations show for A = 0.01 
and 0.0001 that the squared distance is a complicated function of the time. (For A = 1 and 
a > 1 the Qi are larger than one which makes little sense, and for a = and 0.5 at A = 1 
the squared distances are close to t/2; not shown.) We thus look for smaller A at the slopes 
in the log-log plots, i.e. at the effective exponents 

a^s = d\n{Ar{tf)/dt. (4) 

In each case we simulate p = pc, Pc — 0.1, Pc — 0.2, . . . down to pc — 0.5 ^ 0.093. We 
see for short times different slopes in our log-log plots, but for long times the effective 
slopes approach unity: Normal diffusion with mean squared distance proportional to time. 
For a = 2, for an exponentially decaying distribution /(g), and for a WeibuU distribution 
(stretched exponential) the time variations of the effective exponents were similar but less 
pronounced. 

Experimentally more relevant are three instead of two dimensions, and some results are 
shown in FigureOa, rather similar to two dimensions in Figured Now Pc = 0.3116. For very 
small p = 0.0116, squares in Figure [3l we see an overshooting with an effective exponent 
a^s above unity at intermediate times; this is not a statistical fluctuation and shows up 
in all 20 simulated samples (not shown). It also was seen in two dimensions at very low 
p, Figure [3] b. One may call this effect superdiffusion since for more than one order of 
magnitude the exponent is above unity. 

Basically, the positive probability of each barrier to move away and to let through the 
random walker means that for sufficiently long times we always get normal diffusion, a = 1. 
For times which are not long enough to see this moving-away of the barriers, but long 
enough for the walker to explore the whole finite cluster for p < p^ on which it started, 
we have a = 0. For our moderately small A = 0.01 these different regimes cannot be 
reliably separated; that works better for much less mobile barriers: A = 0.0001 in Figure HI 
There the effective exponents are about zero for t ~ 10^, show a maximum but no longer 
overshooting below t ~ 10^, and approach unity above t ~ 10^. 
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IV. DISCUSSION 



In summary, we see a non-monotonic variation of the effective exponents with time, show- 
ing both subdiffusive and super diffusive behavior. Asymptotically, however, the exponent 
always seems to approach unity for t oo. In experiments with more limited variation of 
times, this variation of acs with time could wrongly be interpreted as asymptotic subdiffu- 
sion or asymptotic superdiffusion; long times are needed. 
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